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Abstract. We model energy transfer between two coupled four- level chromophores 
with arbitrarily spaced energy levels. Our analysis takes into account the crucial — 
yet often ignored — process of initial excitation by light that is incident on the 
chromophores. We show that the amount of entanglement generated between the 
chromophores is strongly dependent on the degree of initial excitation as well as the 
inclusion of higher energy levels. We apply our model to the specific example of 
chlorophyll. Our results suggest that an excitation-dependent approach should be 
employed for entanglement studies on multi-level light- harvesting systems even when 
a two-level approximation is valid. 
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1. Introduction 

Photosynthesis is a process utihzed by organisms to produce biomass from raw materials 
such as carbon dioxide and water, using harvested solar energy as the driving force 
(Chapter 12 of [Ij). Central to this process are light-harvesting complexes (LHCs) which 
contain several light-absorbing molecules called chromophores [2l|3l|4]. Photosynthetic 
mechanisms are attracting increasing attention, not only because photosynthesis inspires 
solutions to present renewable-energy needs [51 [6l [71 [HI El [lEl [H] , but also due to recent 
evidence and predictions of non-trivial quantum coherent behaviour in certain LHCs 
[ia[l3l[H[l5l[l6l[I7l[l8l[l9l[2a 

Of additional interest is whether these quantum behaviours play any functional 
role in biology. In non-biological systems, entanglement — a quintessential quantum 
phenomenon — has useful application in enhanced measurements, information processing 
and quantum computation (see Chapter 11 of [25], and [26l |2Il EHl [29l |30l ISI]). It is 
therefore intriguing to ask whether entanglement is present in biological systems and, if 
so, whether organisms utilize entanglement for practical purposes. 

Previous entanglement studies on LHCs have employed two-level models in the 
single-exciton manifold [32l |33l |34l |35j . In this picture, quantum coherence (i.e. nonzero 
off-diagonal elements of the density matrix in the energy eigenbasis) is a necessary 
and sufficient condition for entanglement [36j. An alternative model which treats 
chromophores as quantum harmonic oscillators (systems of an infinite number of equally 
spaced energy levels) excited by a coherent state of light predicts no entanglement in the 
system at all [37] . It is clear that the choice of model, for the LHC and its interaction 
with light, plays a large role in determining the nature of the entanglement generated 
in the system. In this paper, we investigate the features of this role. 

Inspired by known spectroscopy of chlorophylls (Chapter 1 of [38], and |39l [401I4T] ) . 
we consider a model with multiple arbitrarily-spaced energy levels, which we motivate 
in the context of previous models in Section [2} In Section |3| we introduce the formalism 
for our multi-level model. We explicitly take into account a realistic excitation process 
in the preparation of the initial state, rather than assuming that the LHC absorbs 
single-photon pulses. We then evolve the system and calculate the entanglement, 
quantified by the entropy of entanglement, which is generated during energy transfer. 
We present numerical simulations for the entanglement in the system using realistic 
physical parameters in Section |4| and also compare this model with simplified variations 
which consider fewer energy levels, before concluding in Section [5] 

Before we continue, we would like to make a comment on coherence. We understand 
coherence to mean "having a phase relationship". With this definition, it is clear 
that, when used in isolation, the word "coherence" does not mean anything — a phase 
relationship can only exist between two or more properties. Any confusion that 
arises from use of the word is likely to be related to its use in isolation. In this 
paper, we use the word "coherence" in three different contexts. The first is quantum 
coherence in photosynthesis, which refers to coherence (a phase relationship) between 
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energy eigenbasis states of coupled chromophores. The second is in the context of a 
coherent state of hght — a state of indefinite photon-number where the coherence (phase 
relationship) is between the different photon-number states that make up the state. The 
third is in terms of an analogous coherent state of a chromophore — a state with infinitely 
many excited levels — in which the phase relationship is between the energy eigenstates 
of the chromophore. 

2. Modelling light harvesting complexes 

In this section, we summarize previous LHC models and motivate our model. 

2.1. Two-level system in the single- exciton manifold 

By far the most common model for a chromophore is that of a two-level system, 
consisting of an electronic ground state |0) and an excited state |1) |32l [33l [34l [35] . 
Organisms containing these chromophores live at low light conditions under which 
excitation of the chromophore to a higher excited state is unlikely. This is used to 
justify the two-level approximation, which states that the system will only contain at 
most a single excitation; as well as the single-exciton assumption, which states that the 
system will contain only a single excitation. In the single-excitation manifold, the initial 
excited state of an LHC is assumed to have one chromophore in the |1) state and all 
others in |0). These models have predicted long lived entanglement in a chlorophyll- 
containing complex [32] and the Fenna-Mathews-Olson (FMO) complex [33l [34[ [35] . 

However, implicit in the preparation of this initial state is a light-matter interaction 
between a two-level system and a single photon — a highly non-classical state of light. 
This single photon must be absorbed with unit probability to generate a state with 
exactly one excitation, i.e. 

We are at present unaware of any biological mechanism or experiment resulting in 
the absorption of exactly one photon by a chromophore. After interaction with classical 
light like a laser or sunlight, a chromophore would be excited into a superposition or 
statistical mixture of |0) and |1) respectively (Chapter 4 of [25j). This can be thought 
of as allowing for the possibility that light is not absorbed, leaving the system with a 
ground state population. In fact, only around 1% of input sunlight is estimated to be 
absorbed by chlorophyll under natural conditions (see Chapter 5 of [42j). Accordingly, 
the single excitation assumption imposes an excited state that may not correspond to a 
realistic situation. We note that even if it were a single photon that was incident on the 
chromophore, the chromophore would still be excited into a superposition of |0) and |1) 
unless one could guarantee that the photon were absorbed. 

2.2. Quantum harmonic oscillator model 

An alternative model which treats the chromophore as a quantum harmonic oscillator 
(QHO) was introduced by Tiersch et al. [37j. Such a system consists of an infinite 
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number of equally spaced energy levels. Under low-light conditions, the extra levels 
may seem superfluous and the two-level model should suflice, however the rationale for 
the QHO model becomes apparent when considering the light-matter interaction. 

In the dipole and rotating wave approximations, the evolution of the system due 
to the light-matter interaction between a coherent state of light and the QHO is given 
by the displacement operator D{a) = exp{ad'^ — a*a). This immediately gives the 
analogy between the state prepared by this light-matter interaction and the coherent 
state \a) = D{a)\0). An analogy is also made between the FMO complex and a multi- 
armed interferometer. The propagation of a single excitation through the complex is 
analogous to the propagation of a single photon through the interferometer — both will 
clearly lead to entanglement [43j. In contrast, the propagation of the coherent state 
of the chromophore through the complex is analogous to the propagation of a coherent 
state of light through the interferometer, which results in no entanglement being created. 

Tiersch et al show that higher-level truncation of this unentangled state introduces 
the perception of entanglement. A more dramatic increase in the perceived entanglement 
is also observed when the ground state of the coherent state is disregarded, which 
corresponds to working in the single-excitation manifold. 

Although this model realistically considers the process of initial state preparation, 
it only deals with a very special quantum mechanical state and the results may not be 
applicable to systems that deviate from this unique case. In particular, the energy level 
structure of a QHO does not represent that of a chromophore very well. 

2.3. Mutli-level model 

The two models described in the previous sections lead to contradictory results. This 
indicates that initial state preparation and the details of the level structure should 
be taken into account carefully. In our model, we consider a realistic light-matter 
interaction for the state preparation as well as a more complete description of the 
electronic structure of the system itself. 

Certain chromophores do possess higher excitation levels that cannot be neglected 
under typical light conditions. The coherent state of a QHO is a convenient model for 
a many-level system, but it restricts the system to a ladder of equally spaced energy 
levels. The chromophores in question do not have equal energy spacing nor are they 
necessarily in the form of a ladder system. 

A classic example of a multi-level system is chlorophyll in higher plants. The 
absorption spectrum of chlorophyll has two intense bands denoted as Q^ and Soret, 
corresponding to electronic transitions from ground to two separate excited states 
(Chapter 1 of ^3)- Another band denoted as Q^^ is weaker, and thus will not 
be considered in our calculation. Excited state absorption has also been observed 
experimentally in these systems [lH HOI |41] . Another example of multi-level systems 
are carotenoids: chromophores that have at least three observed transitions, i.e. the 
So ^ S2 transition, the ^ Sn transition and the two-photon allowed Sq ^ Si 
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transition (Chapter 1 of [38]). 

With these systems in mind, we develop a more general multi-level model. We 
consider a four-level system with arbitrarily-spaced energy levels, depicted in Figure 
[l] a). In the language of quantum information ^5], this system is referred to as a 
qudit (a rf-level generalization of the quantum bit, or qubit). To model energy transfer, 
we take the simplest case of a coupled dimer consisting of two qudits. The state 
is prepared by interacting a qudit with a pulse of coherent classical light. We note 
that although excitation may not actually occur in the site (chromophore) basis but 
rather in the exciton basis, we will consider the former case in order to put our work in 
the context of previous work on the subject which has considered site-basis excitation 

[aziEaisiEaES]. 




qudit A qudit B 

= -) ^= - 




coupled 



Figure 1. a) Electronic structure of the four- level system (qudit). Transitions from 
|0) to |1) and from |0) to |2) are analogous to transitions leading to the Qy and Soret 
absorption bands in chlorophyl, respectively [38l|39]. The transition from |1) to |3) 
describes excited state absorption (ES A) also witnessed in experiments 011 [40l |4T] . b) 
Two identical uncoupled qudits begin in the ground state, c) qudit A is excited by 



classical light (refer to Section 3.1). d) energy transfer occurs between qudit A and B 
(refer to Section [3^. 



In this paper, we do not consider any decoherence mechanisms which arise from 
the interaction of the chromophores with their environment, nor do we consider line 
broadening or static disorder. These effects are critical for a complete description of 
photosynthetic systems, however, they are beyond the scope of this paper. Our model 
is idealized in this sense and instead, we concentrate on the consequences of different 
excitation processes and energy level structures. 

3. Theory 

To model the entanglement generated during energy transfer in a multi-chromophoric 
system, we study the simplest case of two coupled qudits. We prepare the initial state 
by considering the light-matter interaction between classical light and one qudit, and 
then calculate the evolution of the entanglement between the two coupled qudits. 
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3.1. Initial excitation 



For excitation of chromophore A in the site basis, we first assume that the chromophores 
are uncoupled. We take the total Hamiltonian for the qudit-light system to be 

Hi = Ha + , (1) 

where is the Hamiltonian of the qudit given by 



= ^hwn\n){n\ 

n=0 



(2) 



where /kJo is the ground state energy and /icji_3 are energies of the three excited states. 

As the size of a chlorophyll molecule is roughly Inm — much shorter than the 
wavelength of visible light — the dipole approximation for absorption can be invoked 
(Chapter 3.6 of [46j), giving the light-matter interaction Hamiltonian 

^Im = -E{t) . I (3) 

where d is the dipole moment operator and E = E^^\t) + E^~\t) is a classical electric 
field given by 

E^^\t) = Eoe'^'^'e, (4) 

where £"0 is the electric field amplitude, cjl is the carrier frequency and e is the 
polarization vector. We can define a co-rotating frame, where the transformation is 
given by the unitary operator U such that 

\<Pcr) =U\<P), (5) 

where 10) is the state in the Schrodinger picture. We aim to find a time-independent 
Hamiltonian such that the dynamics of the system is given by (Chapter 9 of p7] ) 

\Mt)) =e-t^-*|</>e.(0)), (6) 

where 



Her = U^Hi{t)U - i^U^^U . 



(7) 



We find that the operator 



u 



V 
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-itUJL 



\ 



J 



(8) 



gives the following time-independent Hamiltonian in the rotating wave approximation 



Her = 



v 





Eodoi 


EodQ2 





EodiQ 


huji 





Eodis 


Eod2o 





huj2 








Eodsi 








(9) 



where dnm = ^' {n\d\m) and we have made use of the experimental observation that in 
chlorophyll, the most physically relevant transitions are those illustrated in Figure [l] a). 
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3.2. Inter- qudit excitation energy transfer (EET) 

To model energy transfer, we consider two coupled qudits. We take the Hamiltonian for 
this system to be 

H2 = Hm + He , (10) 

where = -^d ® Ib + Ia ® -^d is the Hamiltonian for the free qudits, where is 
defined in ([2]), and 

He = (jio,oi|10)(01| + J2o,o2|20)(02| + Ji3,3i|13) (31| 

+ J20,oi|20)(01| + Jio,o2|10)(02| + Ji2,3o|12)(30| (11) 
+ Jii,3o|ll)(30| + Jn,o3|ll)(03|) +h.c. 



This coupling Hamiltonian is constructed according to the observation that only 
transitions between |0) and |1); |0) and |2); and |1) and |3) can occur in the systems 
under consideration [38l [40l El |44] . 

3.3. The final state 

To calculate the final state, we begin with both uncoupled chromophores in the ground 
state |'0(O))a,b = |0)a|0)b as indicated in Figure [l]b). After interaction of qudit A with 
light, as indicated in Figure [l]c), the state is 

|^in(r))a = f/te-^^-^[/|0)a, (12) 

where U is defined in ([s]) and Her is given in We take |'0in(r))a,b = |'0in(r))a|O)b 
to be the initial state in terms of the impending energy transfer between the coupled 
qudits. We now let the qudits couple to each other to allow for the energy transfer, as 
indicated in Figure [l]d). The final state is 

|^(t))a,b = e-t^^*|^in(r))a,b, (13) 

where t ^ t + T since t ^ T, T is the duration of the light-matter interaction, and H2 



is defined in (10). 



3.4' Measure of entanglement 
Any pure bipartite state 

3 3 

\Mt)) = 5^5Zc^n|m)|n), (14) 

can be decomposed into what is known as the Schmidt-decomposed form 

3 

k=0 

where {lo^j^^)} and {|/3/c)} form orthonormal bases and the diagonal matrix s is found by 
performing a singular value decomposition of the matrix c (Section 20.2 of [48j). 
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We use the entropy of entanglement, defined as the von Neumann entropy of one of 
the reduced states of a bipartite system, as the measure of the entanglement over time 
for the two-qudit system. In terms of the Schmidt values this is ^48j i49j 

3 

Emm =-^4^og,si (16) 

k=0 

The entropy of entanglement ranges from zero for a product state to log2 for a 
maximally entangled state of two A^-state particles, which in our case is log2 4 = 2. 

3.5. Parameters 

We estimate reasonable physical parameters by considering an analogous photosynthetic 
system of two weakly-coupled pairs of strongly-coupled qubits [50j and using numerical 



values from [5l]. For a detailed calculation, refer to Appendix A Our estimated 
parameters are summarized in Table [TJ 

Level Energy Dipole Strength Coupling 
^0 = 

uji=u) \dio\ 

UJ2 = l.OAuj \d2o\ 



'31 



= d 


Jwfii 


= J 


= OMd 


^20,02 


= 0.50J 




■^20,01 


= -0.67J 




^10,02 


= 0.72J 


= d 


■^13,31 


= 0.90J * 




^11,30 


= ^11,03 = 




■^12,30 


= 0.76J * 



0.81J 



Table 1. Estimated numerical values for relative energy spacings and couplings. * 
denotes values that are empirically assigned by assuming that J is slightly larger than 
-^13,317 which is slightly larger than Jii,3o = -^11,03, which is in turn larger than J12, so- 
All other values are estimated according to [50], [51] and [12 . Typical values are 
uj = 2.99 X 10^^s-\ d = 5D and J = -1.25 x lO^^J (see [Appendix A[). 



To facilitate interpretation of the numerical results, we introduce four convenient 
dimensionless parameters, given in Table [2j 
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Parameter 



Interpretation 



7 = E^dT/h 

8 (jjT j(juh/Eod 

72 = Jt/h 

82 = (jJt = 72Cl;/i/J 



degree of initial excitation 



coupling strength and/or evolution time 



Table 2. Summary of dimensionless parameters. £^0 is the electric field amplitude, d 
is related to the dipole moments (see Table IT]), T is the light-matter interaction time, 
u is related to the level energies (see Table Q, J is related to the coupling strengths 
(see Table [T]) and t is the evolution time, where t ^ T. We study how entanglement 
changes as 7 and 72 vary while 5 is kept fixed. Note that 82 varies as a function of 72 



as 82 = 72i^/J. A typical value is 7 = 0.41 (see Appendix B). 



4. Results and discussion 

In this section, we numerically investigate the entanglement between the two 
chromophores as a function of the degree of excitation and the dynamics of the state. 

4.1. Excitation dependence 

To investigate the excitation dependence on the prepared initial state, we plot 
probabilities = |('0in(r)|a|^)aP as a function of 7 in Figure [2| The maximum 
entanglement between the two chromophores E'max = niax(£'[|'0(t))a,b]) as a function 
of 7 is also shown in Figure [2} We find that E'max is strongly dependent on the degree 
of the initial excitation. 




no truncation 

single-exciton 

manifold 
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Figure 2. Probability Pn = |('0in(^)|a|^)aP foi" qudit A and the resultant maximum 
entanglement max(£^[|?/^(t))a,b]) after inter-qudit coupling. Note that, in accordance 
with the definitions in Table |2j a more intense light field leads to a larger 7. Within 
< 7 < 1 which is close to the typical experimental condition of 7 = 0.41, increasing 
7 gives a higher probability to be in the excited states and subsequently a larger 
maximum entanglement between qudits. 
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For small values of 7, e.g. those corresponding to typical light conditions such as 
7 ^ 0.41 (as calculated in Appendix B), E'max increases with the level of initial excitation. 
Comparison of this with the maximum entropy of entanglement when assuming a single- 
exciton manifold, which is equivalent to assuming an initial state of |'0in(r))ab = |10)ab5 
is also shown in Figure [2] 

Incidentally, at 7 = 3, qudit A has a probability of pi ^ 0.9. When qudit A couples 
to qudit B, the |2) and |3) contributions in either qudit remain low and entanglement 
can reach E'max ^ 1, also shown in Figure [2j This special regime resembles the case of 
the single-excitation manifold model. In the other limit when 7^0, the entanglement 
vanishes, just as it does in the QHO model [37j. These distinct entanglement outcomes 
depict the dependence of entanglement on light-matter interaction. 

Figure [3] illustrates the dramatic dependence of the entropy of entanglement on 72 
as well as 7. For small values of 7 the entropy of entanglement increases with the level 
of initial excitation for all values of 72. 




Figure 3. Dependence of the entropy of entanglement on the degree of initial 
excitation. In the regime of < 7 < 1, a more highly-excited qudit generates a 
larger entropy of entanglement when it is coupled to another ground-state qudit. The 
cross section at 7 = 0.41 is plotted in Figure |4] in cyan. Clearly, the degree of initial 
excitation quantified by 7 largely determines the entanglement profile, and should thus 
be considered when examining entanglement in LHCs. 



4.2. Effect of neglecting higher levels 

To evaluate the effectiveness of simplified versions of our model, we consider a number of 
truncations of the four-level system. The Hamiltonians have been modified accordingly 
to incorporate the simplification. The entropy of entanglement as a function of 72, for 
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no truncation, 7 = 0.41 
no truncation, 7 = 3 
truncate |3), 7 = 0.41 
truncate |2) and |3), 7 = 0.41 
single-exciton manifold 



Figure 4. The entropy of entanglement for the four-level system for two different 
values of 7; the three-level and two-level simplifications; as well as the case of the 
single-exciton manifold. 

each truncated system compared with the untruncated case, is shown in Figure |4} The 
case of the single-excitation manifold is also shown. Neglecting the highest excited state 
(i.e. the |3) state) does not significantly change the entanglement profile compared to 
the four-level case, as can be expected from the fact that the |3) state already has 
minimal contribution even before being neglected. Accordingly, we consider the three- 
level simplification to be a good approximation. 

Subsequent neglection of the |2) level reduces the maximum perceived entanglement 
by approximately half. A two-level simplification may therefore provide a lower 
bound for the estimate of entanglement for a pair of four-level systems with a V- 
like configuration. This is in contrast with the QHO model (a ladder system) where 
neglecting higher levels always increases the perceived entanglement [37], emphasizing 
the importance of careful consideration of the energy level structure of the system. 
Notice that even in the case of a two-level approximation that considers the light-matter 
interaction, the entanglement profile is dramatically different to that of the single-exciton 
manifold, showing that one does not imply the other. 

It is perhaps surprising that neglecting the possibility of Soret absorption has a 
much more dramatic effect on the entanglement than does neglecting the possibility of 
excited-state absorption. 

5. Conclusion 

We provide a general model for a pair of four-level systems in a V-like configuration 
with arbitrary energy spacings, and input numerical parameters to obtain physically 
reasonable values for the entropy of entanglement. 

We have demonstrated that the entanglement profile largely depends on the degree 
of initial excitation, which incorporates the light intensity, field duration and transition 
dipoles. Within a realistic range of physical parameters, a more excited chromophore 
generally produces larger entanglement upon interaction with the other initially un- 
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excited chromophore. Under extreme degrees of initial excitation, our model can 
produce a maximum entanglement either as high as in a two-level model in the single- 
exciton manifold or as low as in a harmonic oscillator model. This dramatic difference 
demonstrates the importance of considering the light-matter interaction. 

Surprisingly, if a system has four levels in this particular configuration, as in 
chlorophylls, the two-level simplification may provide a lower estimate of entanglement. 
Nevertheless, the excitation process should be considered to determine the initial state. 
Otherwise, the single-exciton manifold may significantly overestimate the entanglement. 

We did not consider any decoherence mechanisms which arise from the interaction 
of chromophores with their environment, nor did we consider line broadening or static 
disorder. It will be interesting to incorporate our findings into a decoherence model to 
study the entanglement in a photosynthetic complex like chlorophyll — one of the most 
common yet sophisticated systems in nature. 
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Appendix A. Parameter estimation 



For chlorophylls, numerical values of certain couplings, such as ^20,02 in (12), are not 
available. Thus we seek estimates of parameters that are within the physical order 
of magnitude for photosynthetic systems. It has been shown in [50j that a system 
consisting of four qubits with the following properties — qubits Ai and A2 are strongly 
coupled, qubits Bi and B2 are strongly coupled, while qubits A^ are weakly coupled to 
qubits B^ — can be treated as a two-qudit system. 

These conditions are satisfied by the photosynthetic protein Phycocyanin 645, where 
we define the chromophores MBV 19 A, PCB 158 C, MBV 19 B, and PCB 158 D as 
qubits Ai, A2, Bi and B2, respectively (see Table S.l in [52]). Chromophores are often 
referred to as sites in the literature. These four sites effectively form a two-qudit system: 
sites Ai and A2 couple to form qudit A while sites Bi and B2 couple to form qudit B. 
We split the Frenkel exciton Hamiltonian (which is given in the single-exciton subspace) 
for the combined four-site system into two parts: 




J^AB= " ^ \ +V, (A.l) 



V incorporates the weak couplings between Ai and 5^, and is given in (A. 16). Ha and 
Hb are the Frenkel exciton Hamiltonians for the combined Ai and A2 system and the 
combined Bi and B2 system respectively: 

Ha = I 1 , (A.2) 



VA2A1 



2 
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Diagonalizing Ha and Hb gives eigenvectors 

= «ii|eAi) + ai2|eA2> , 
[2) A = a2i|eAi) + a22|eA2) , 



(A.3) 



(A.4) 
(A.5) 



and 



= ^iile^i) + fei2|eB2) , (A. 6) 

|2)B = &2i|eBi)+&22|eB2), (A.7) 

with corresponding eigenvalues cj^,!, ^a,2 and cj^^i, cc;^,2 respectively, i.e. the energies 
of the two lower excited states of each qudit. 

Energies of the third excited states |3)a and \?>)b are estimated to be uja,"^ = "^^a.i 
and (jJb^?, = Scj^^i based on a typical chlorophyll molecule where the energy of the 
|1) 1 3) transition (excited state absorption) is close to the energy of the |0) |1) 
transition (Q^; absorption) [44|, HOI [41] . For the ground states |0)a and |0)b, we define 



Starting from the values 
Ha = 



16050 
-87 



and 



16373 
86 



-87 
15808 



15889 



taken from [52], we find that wa,i ~ wb,i and wa,2 
= cja,i = uJb,i = 2.99 X lO^^s"^ , 

U2 = L0a,2 = 1^5,2 = 3.11 X lO^^S"^ . 

Couplings between the qudits are given by 



-^10,01 


{MaV\1)b 


-^20,02 


{2\aV\2)b 


-^10,02 


{1\aV\2)b 


-^20,01 


{2\aV\1)b 



{A.i 



(A.9) 

Wb,2' Thus for simplicity we take 

(A.IO) 
(A.ll) 

(A.12) 
(A.13) 
(A.14) 
(A.15) 



where the states are defined in (A.4 A.7) and 

/ 4 -3 \ 

3 8 

4 3 

-3 8 



V 



\ 



cm 



is based on parameters given in [52j. This gives J = Jio,oi 
couplings are given as ratios in Table [T] 



(A.16) 



-1.25 X 10^2 J and other 
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The transition dipoles for qudit A are given by [5^ 



dB,10 
dB,20 



Ad\0)A. 

Ad\0)A, 

Bd\0)B , 



where the states are defined in ( |A.4 A.7) and 

(eAi|rf|0) = (-1.42,4.54,-13.70) D, 
(eAM^) = (13.58,3.53,1.78) D, 
(eBi|rf|0) = (1.50,2.60,-14.20) D, 
{cbMI^) = (4.98,-12.51,-3.81) D 

are the dipole moments for qubit Ai, A2, Bi and B2 given in [51]. We take \djk\ = 
\dAjk\ = \dB,jk\^ where A; = — 3. We also take |rf3i| = |rfio| based on available data 
which show that the transition strength of the excited state absorption is close to that 
of the absorption [44[ HOI [41]. These physical parameters give \dio\ = 5D and other 
dipole moments are given as ratios in Table [T] 



(A.17) 
(A.18) 
(A.19) 
(A.20) 

(A.21) 
(A.22) 
(A.23) 
(A.24) 



Appendix B. Light intensity calculation 

The electric field Eq in a laser pulse can be estimated by 

2 V ' 

where eo and is the permittivity of free space, W is the energy of the pulse and V = cTA 
is the beam volume given by the speed of light c, pulse duration T and beam cross-section 
A. From this, we can define a dimensionless parameter 



(B.l) 



^ h h\ cAe^ ^ ^ 

in terms of known parameters. Pulses used to illuminate photosynthetic systems have 
typical values oiW = 5nJ, T = lOfs and A = 25007r/xm^. We also use = 5D (calculated 



in Appendix A) to estimate a typical value of 7 = 0.41. 
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